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Abstract 

This paper is devoted to the analysis of some uniqueness properties 
of a classical reaction-diffusion equation of Fisher-KPP type, coming 
from population dynamics in heterogeneous environments. We work 
in a one-dimensional interval (a, b) and we assume a nonlinear term of 
the form u {n{x) — ju) where fj, belongs to a fixed subset of C°([a, b]). 
We prove that the knowledge of u at t = and of m, at a single 
point xq and for small times t € (0, e) is sufiicient to completely de- 
termine the couple (u(t,x), iJ,{x)) provided 7 is known. Additionally, 
if Uxxit,XQ) is also measured for t G (0,6), the triplet x), /i(a;), 7) 
is also completely determined. Those analytical results are completed 
with numerical simulations which show that, in practice, measurements 
of u and Ux at a single point xo (and for t S (0,£)) are sufficient to 
obtain a good approximation of the coefficient ^i{x). These numerical 
simulations also show that the measurement of the derivative Ux is 
essential in order to accurately determine n{x). 

Keywords: reaction-diffusion • heterogeneous media • uniqueness • inverse 
problem 



1 Introduction and ecological background 

Reaction-diffusion models (hereafter RD models), although they sometimes 
bear on simplistic assumptions such as infinite velocity assumption and com- 
pletely random motion of animals [Tj, are not in disagreement with certain 



1 



dispersal properties of populations observed in natural as well as experimen- 
tal ecological systems, at least qualitatively [11 [3l IH [5]. In fact, since the 
work of Skellam [6], RD theory has been the main analytical framework to 
study spatial spread of biological organisms, partly because it benefits from 
a well-developed mathematical theory. 

The idea of modeling population dynamics with such models has emerged 

th 

at the beginning of the 20 century, with random walk theories of organ- 
isms, introduced by Pearson and Blakeman [7]. Then, Fisher [B] and Kol- 
mogorov, Petrovsky, Piskunov [9] independently used a reaction-diffusion 
equation as a model for population genetics. The corresponding equation is 

du d'^ u 

- D-^ = u (/i - 7n), t>0, x£{a,b) C R, (1.1) 

where u = u{t, x) is the population density at time t and space position x, 
D is the diffusion coefficient, and fj, and 7 respectively correspond to the 
constant intrinsic growth rate and intraspecific competition coefficients. In 
the 80's, this model has been extended to heterogeneous environments by 
Shigesada et al. [TU]. The corresponding model is of the type: 

The coefhcients /-f(x) and 7(x) now depend on the space variable x and can 
therefore include some effects of environmental heterogeneity. More recently, 
this model revealed that the heterogeneous character of the environment 
played an essential role on species persistence and spreading, in the sense 
that for different spatial configurations of the environment, a population 
can survive or become extinct and spread at different speeds, depending 
on the habitat spatial structure ([2], [11], [IS], [13], [14] ,[I5], [16]). Thus, 
determining the coefficients in model ()1.2p is an important question, even 
for areas other than ecology (see jl7j and references therein). 

In this paper, we focus on the case of constant coefficients D and 7: 

du d'^u 

'Ol ~ ^ ^ ^^^^^ -lu), t > 0, X e (a, b), (1.3) 

and we address the question of the uniqueness of couples (n, /i(x)) and triples 
(«, /Li(x), 7) satisfying ()1.2p . given a localized measurement of u. 

Uniqueness results of this type have been obtained for reaction-diffusion 
models, through the Lipschtiz stability of the coefficient with respect to the 
solution u. Lipschtiz stability is generally obtained by using the method of 
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Carleman estimates [18]. Several publications starting from the paper by 
Isakov [19] and including the recent overview of the method of Carleman 
estimates applied to inverse coefficients problems [20] provide results for the 
case of multiple measurements. The particular problem of the uniqueness of 
the couple {u, fi{x)) satisfying (|1.3p given such multiple measurements has 
been investigated, together with Lipschtiz stability, in a previous work [21] . 
Placing ourselves in a bounded domain Q of M.^ with Dirichlet boundary 
conditions, we had to use the following measurements: (i) the density u{0, x) 
in at t = 0; (ii) the density u{t,x) for (t,x) € (toj^i) ^ for some times 
< to < ^1 and a subset uj CC ^2; (iii) the density u{6,x) for all x € i}, at 
some time 6 € {tQ,ti). 

Although the result of [21] allows to determine using partial mea- 
surements of u{t, x), assumption (iii) implies that u has to be known in the 
whole set fi. This last measurement (iii) is a key assumption in several other 
papers on uniqueness and stability of solutions to parabolic equations with 
respect to parameters (see Imanuvilov and Yamamoto [22], Yamamoto and 
Zou [23], Belassoued and Yamamoto [23] for scalar equations and Cristofol, 
Gaitan and Ramoul [25] or Benabdallah, Cristofol, Gaitan and Yamamoto 
[26] for systems). 

Here, contrarily to previous results obtained for this type of reaction- 
diffusion models, there are some regions in (a, b) where u is never measured: 
we only require to know (i') the density u{0, x) in (a, b) at t = and (ii') the 

du 

density n(t,xo) and its spatial derivative — (t, xq) for t S (0,e) and some 

point xq in (a, b) (see Remark 12.51 for a particular example of hypothesis 
(ii')). Thus a measurement of type (iii) is no more necessary. Furthermore, 
we show simultaneous uniqueness of two coefficients n{x) and 7 provided 

that measurements of the second derivative Tr-^itjXo) are available. 

OX'^ 

Our paper is organized as follows: in the next section, we give precise 
statements of our hypotheses and results; Section [3] is then dedicated to the 
proof of the results. Section [3] is devoted to the description of numerical 
examples illustrating how the coefficient /u(x) can be approached using mea- 
sures of the type (i') and (ii'). Those results are further discussed in Section 

m 
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2 Hypotheses and main results 

Let (a, b) be an interval in M. We consider the problem: 

' f -Z)0 = n(M(x)-7n), t>0, x G {a,b), 
^ aiu{t,a)-/3i^{t,a)=0,t>0, , . 

«2u(t,5) + /32f =0, t >0, ^ 
^ m(0, x) = Ui{x), X € (a, 6). 

Our hypotheses on the coefficients are the following. Firstly, we assume 
that: 

/i € M := G C^'^{[a,b]) such that ^ is piecewise analytic on (a, 6)}, 

(2.4) 

for some 7] € (0, 1]. The space C^'"^ corresponds to Holder continuous func- 
tions with exponent r] (see e.g. [2Z])- A function ^p G C^''^{[a,b]) is called 
piecewise analytic if it exists n > and an increasing sequence {ik)i<k<n 
such that ii = a, in = b, and 

n-l 

for all X G (a, 6), ip{x) = ^X[i,,i,+i){x)^Pj{x), 

i=i 

for some analytic functions <pj, defined on the intervals and where 

X[ij,ij+i) the characteristic functions of the intervals for j = 

l,...,n-l. 

We also assume that 7 is a positive constant and that the boundary 
coefficients satisfy: 

ai,a2,Pi,f^2 > with oi + /3i > and a2 + > 0. (2.5) 
We furthermore make the following hypotheses on the initial condition: 

Ui>0, Ui^O and u^ G C^'^^iia, b]), (2.6) 

for some r] in (0, 1), that is Ui is a function such that u'- is Holder contin- 
uous. In addition to that, we assume the following compatibility conditions: 

aiUi{a) - /3iu-(a) = 0, a2Ui{b) + (i2u'i{b) = 0, Sfs^u'-^a) = 0, 8p^u'l{b) = 0, 

(2.7) 

where 5y is verifies: 5q = 1 and 5^ = if y 7^ 0. We also need to assume 
that: 

measure({x G (a, 6), Ui{x) = 0}) = 0. (2-8) 
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Under the assumptions ()2.4p - (j2.7p . for each ^ ^ M and 7 > 0, the 
problem {Vf_i,'y) has a unique solution u € C^'^^2([0) +00) x [c^)^]) the 
derivatives up to order two in x and order one in t are Holder continuous, 
see [271 EH] for a definition of Holder continuity). Existence, uniqueness and 
regularity of the solution u are classical. See e.g. [251 Ch. 1]. 

Let us state our main results: 

Theorem 2.1. Let fJ-,flG M, and 7 > 0, and assume that the solutions u 
and u to ('^^,7) and (7^/2,7) satisfy, at some xq € (a, 6), and for some e > 
and all t in (0, e).- 

u{t,xo) = u{t,xo), (2.9) 
^(i,^o) = ^(*>^o). (2.10) 

Assume furthermore that 

Uiix^) ^0 or Q^it, xo) = -Q^{t^ Xq) for t £ (0, e). (2.11) 

Then, we have fj, = fi on [a,b] and consequently u = u in [0, +00) x [a,b]. 
If Pi > (resp. (32 > 0), this statement remains true when xq = a (resp. 
Xq = b). 

Remark 2.2. This result remains valid if j = 7(2;) is a given, positive 
function in C^''^{[a,b]). 

However, the conclusion of Theorem 12.11 is not true in general without 
the assumption (|2.10p : 

Proposition 2.3. Let fi M and 7 > 0. Assume that a\ = 02 and /3i = (^2 

and that Ui is symmetric with respect to x = {a + b)/2. Let fx := ^{b—{x — a)) 
for X G [a,b]. Then, the solutions u and u to (^^,7) and (Vp,,^) satisfy 
u{t, ^) = u{t, ^) for allt>0. 

Under an additional assumption on the initial condition m, we are able 
to obtain a uniqueness result for triples (u, /x, 7): 

Theorem 2.4. Let G M, and 7,7 > 0. Assume that, at some xq € 
{a,b), Ui{xo) = 0. Assume furthermore that the solutions u and u to (Vfj,,^) 
and (Vfi^^) satisfy, for some e > and for all t in (0, e); 

u{t,xo) = il{t,xo), (2.12) 

du , . du , . , . 

^(i.^o) = ^(i'^o), (2.13) 

9^(*'^o) = Q^it,xo). (2.14) 
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Then, we have fi = fl on [a, b] and 7 = 7. Consequently u = u in [0, +00) x 
[a,b]. If Pi > Q (resp. (^2 > this statement remains true for xq = a 
(resp. Xq = b). 

Remarks 2.5. • A particular example where hypotheses \2.9\2.lJ\j of 
Theorem \2. 1\ ( resp. hypotheses \2.12^^2A^ of Theorem \2.4^ are fulfilled 
is whenever, for some subset oj of (a, b), u{t, x) = u{t, x) for t G (0, e) 
and all X LJ (resp. xq € a; and u{t,x) = u{t,x) in (0,e) x ojj. 
Note that, under this hypothesis, the previous results 121^ did not imply 
uniqueness; indeed, an additional assumption of type (Hi) was required 
(cf. the introduction section). 

• The uniqueness result of Theorem \2.4\ cannot be adapted to the sta- 
tionary equation associated to (V^^-y): —p" = p{fj.{x) — ^p) (see e.g. 

for the existence and uniqueness of the stationary state p > 0). 
Indeed, for any r € (0, 1), setting ji = ^ — rjp and 7 = (1 — t)j, 
we obtain —p" = p{jl{x) — 'yp), whereas fi ^ n and 7^7. Thus, a 
measurement of p, even on the whole interval [a,b], does not provide 
a unique couple (/U,7). 

• The subset M 0/ C'^'^([a, 6]) made of piecewise analytic functions is 
much larger than the set of analytic functions on [a,b]. It indeed con- 
tains some functions whose regularity is not higher than C^'"^ , and 
some functions which are constant on some subsets of [a, b] . Our re- 
sults hold true if M is replaced by any subset M' of C^'^{[a,b]) such 
that for any couple of elements in M' , the subset of [a, b] where these 
two elements intersect has a finite number of connected components. 



Let fi,fl € M, and 7,7 > 0. Let u be the solution to {V^i^-f) and u the 
sokition to {Vjx^^). We set 



3 Proofs 



U := u — u and m := fi — jl. 



The function U satisfies: 




flU — jU {u + u) + u{m — u{'~f — 7)) 



(3.15) 



for t > and x G (a, b), and 



aiU{t,a)-/3i^{t,a) = 0, a2U{t,b) + /32^{t,b) = 0, t>0 
U{0,x) = 0, xe {a,b). 



(3.16) 
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Proof of Theorem 12.11 In that case 7 = 7- Equation ()3.15p then reduces 

to 

dU d'^U 

- D-^ = flU - -fU{u + u) + um. (3.17) 
Step 1: We prove that m{xQ) = 0. 

It foUows from hypothesis (j2.9p that, for ah t G [0,e), U{t,xo) = and 
thereby, 

dU 

— {t,xo) = for all t G [0,e). 

If lij(xo) 7^ 0, then, since C/(0, •) = we deduce from p.l7p applied at t = 
and X = xq that Ui{xo)m{xo) = 0, and therefore m{xQ) = 0. 

If Ui{xo) = 0, from (IZTTT) . we have ^{t,xo) = for aU t € [0,e). 
Applying equation (13.17P at t = e/2 and x = xq, we get 



u (I'^o) m{xo) = 0. 



If xq E {a,b), the strong parabolic maximum principle (Corollary 15. 2p ap- 
plied to u implies that u{e/2,xo) > 0. As a consequence we again get 
m(xo) = 0. Lastly, if xq = a and /3i > the Hopf's Lemma applied to 
u again implies that u{£/2,xq) > 0. Indeed, assume on the contrary that 
u{£/2,xq) = u{e/2,a) =0. The boundary condition aiu(t, a) — /3i|^(t, a) = 
implies: 

^ du f£ \ 

=0, 



dx V2' 

which is impossible from Hopf's Lemma (Corollary 15.21 and Theorem 15. II (b) 
and (c)). Thus 'u(e/2,xo) > and, again, m(xo) = 0. A similar argument 
holds for Xq = b, whenever (^2 > 0. 

Under the assumptions of Theorem l2.H we therefore always obtain m(xo) ■ 

0. 

Step 2: We prove that m = 0. 
Let us now set 

hi := sup{x € [xo,6] s.t. m has a constant sign on [xo,x]}. 

By "constant sign" we mean that either m > on [xq, x] or m < on [xq, x]. 
Then, four possibilities may arise: 
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• (i) m = on [xq, bi] and 61 < b, 

• (ii) ?Ti > on [xo,foi] and it exists xi G (2:0,61) such that m(xi) > 0, 

• (iii) m < on [xo,6i] and it exists xi S (xo,fei) such that m{xi) < 0, 

• (iv) 61 = b, and m = on [xq, b]. 

Assume (i). Then, by definition of 61, there exists a decreasing sequence 
Uk ^ bi, yk > bi, such that \m{yk)\ > for all k > 0. Assume that it 
exists ko such that |m(x)| > for all x € (6i,2/fco)- continuity, m does 
not change sign in {bi^yko)-, and therefore in [xcy^o]- This contradicts the 
definition of 61. Thus, 

for all /c, it exists S (6i,yfe) such that m{zk) = 0. (3.18) 

Since /U and fi belong to M, the function m also belongs to M and is therefore 
piecewise analytic on (a, 6). Thus, the set {x G (a, 6) s.t. m{x) = 0} has a 
finite number of connected components. This contradicts ()3.18p and rules 
out possibility (i). 

Now assume (ii). By continuity of m, and from hypothesis (12. 8p on Ui, 
we can assume that Ui{xi) > 0. Since m(xi) > and U{0, •) = 0, it follows 
from (IHTT]) that 

dU , , / x / \ 
— (0,xi) = Ui{xx)m[x\) > 0. 

Thus, for ei > small enough, U{t, xi) > for i € (0, ei]. As a consequence, 
U satisfies: 

^ ^-D^-{il-ju--fu)U>0, ie(0,ei], xeixo,xi), 

U{t, xo) = and U{t, xi) > 0, t e (0, ei], (3.19) 
^ U{0,x) = 0, X G (xo,xi). 

Moreover, 

Lemma 3.1. We have U{t,x) > in (0,ei) x (xo,xi). 

Proof of Lemma Q Set W = f/e"^*, for some A > large enough 
such that c{t, x) := fi — — ju — X < in (xq, xi). The function W satisfies 

dW d'^W 

— - D—^ - c{t,x)W > 0, t G (0,ei], x G (xo,xi). 
at ax^ 

Assume that it exists a point (t*, x*) in (0,ei) x (xq, xi) such that U{t*,x*) < 
0. Then, since W{t,xo) = and W(t,xi) > for t G (0,ei), and since 
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W{0,x) = U{0,x) = 0, W admits a minimum m* < in (0,ei] x (xq,xi). 
Theorem 15. II (a) applied to W implies that W = m* < on [0, ei] x [xq, xi], 
which is impossible. Thus W > in [0, ei] x [xq, xi]. Theorem 15. II (a) and (c) 
then implies that W > and consequently U{t,x) > in (0,ei) x (xq,xi). 
□ 

Since U{t,XQ) = 0, the Hopf's lemma (Theorem 15.11 (b) and (c)) also 
dU 

implies that -—(t,xo) > for all t G (0, ei). This contradicts hypothesis 
ox 

(j2.10p . Possibility (ii) can therefore be ruled out. 

Applying the same arguments to —U, possibility (iii) can also be rejected. 
Finally, only (iv) remains. 

Setting 

ai := inf{x € [a,xo] s.t. m has a constant sign on [xjXq]}; 

the same argument as above shows that ai = a and m = on [a, xq]. Thus, 
finally, m = on [a, b] and this concludes the proof of Theorem 12.11 □ 

Proof of Theorem 12.41 From the assumptions ()2.12p and ()2.14p of Theo- 
rem 12.41 equation ()3.15p at x = xq reduces to 

u{t, xq) (m(xo) - u{t, xo){'y - 7)) = for t G [0, e). 

If Xq G (a, b), the strong parabolic maximum principle (Corollarv l5.2p implies 
that u{t,xo) > for all t > 0. This remains true if xq = a (if /3i > 0) or 
Xq = b (if (^2 > 0); cf. the proof of Theorem 12.11 We therefore get: 

m{xo) = u{t, xo)(7 - 7) for t G (0, e). (3.20) 

From the continuity of t u{t, xq) up to t = 0, we have m{xo) = iii(xo)(7 — 
7). Thus, Ui{xo) = implies that m(xo) = which in turns implies from 
()3.20p . and since u{t, xq) > for t > 0, that 7 = 7. The end of the proof is 
therefore similar to that of Theorem 12.11 □ 

Remark 3.2. Extension of the arguments used in the previous proof to 
higher dimensions is not straightforward. Indeed, placing ourselves in a 
bounded domain Q. o/M^, with N > 2, we may consider the largest region 
Qi in il., containing xq and such that m has a constant sign in Oi. Consider 
in the above proof the possibility (HJn (instead of (ii)): m > on Oi and it 
exists xi £ Q such that m{xi) > 0. Then it exists a subset uJi ofQi, such 
that Xq € duJi and u(t,x) > on a portion of duji. However, we cannot 
assert that U{t,x) > on dooi, and (ii)N can therefore not be ruled out as 
we did for (ii). 
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Proof of Proposition 12. 3t Under the assumptions of Proposition 12.31 we 
observe that u{t, b — (x — a)) is a solution of {V^^^). In particular, by unique- 
ness, we have 

u{t, x) = u{t, b — (x — a)), for all t > and x G [a, b]. 

It follows that u{t, ^) = u{t, ^) for all t > 0. □ 



4 Numerical computations 

The purpose of this section is to verify numerically that the measurements 
(j2. 9112. 10"]) of Theorem 12.11 allow to obtain a good approximation of the co- 
efficient fi{x), when 7 is known. 

Assuming that /i belongs to a finite-dimensional subspace E C M and 
measuring the distance between the measurements of the solutions of (V^^^) 
and {Vfi^-y) through the function 

Gf,{fi) = \\u{-,xo) - n(-,3;o)||L2(o,e) + ll^(-'^o) - -g^{-,Xo)\\L2{0,e), 

we look for the coefficient fJ,{x) as a minimizer of the function G^. Indeed, 
Gfj^{fi) = and, from Theorem 12. 1|, this is the unique global minimum of 
in M. 

Solving (V^^j) by a numerical method (see Appendix B) gives an ap- 
proximate solution u^. In our numerical tests, we therefore replace by 
the discretized functional 

Quj^ Qu}^ 

G^{jj) := ||n''(-,xo) - u^{-,xo)\\L^o,e) + ~ -g^i-,Xo)\\L2{0,s)- 

Remark 4.1. Since (V^^j) and (Vfi^-y) are solved with the same (determin- 
istic) numerical method, we have G^{^) = 0. Thus fi is a global minimizer 
of G^. However, this minimizer might not be unique. 

4.1 State space E 

We fix (a, b) = (0, 1) and we assume that the function /i belongs to a subspace 
E C M defined by: 

E:=h€ G^'^iiO, 1]) I 3 {h,)o<^<n e M"+i, ^(x) = ^hi-j {{n - 2) {x - q)) on [0, 1] 

I i=0 

i-1 J V N f expf^^V ifxG(-2,2), 
with Ci = ^ and ^(2;) = <^ ^ \x^-^J ' v > 

otherwise. 
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4.2 Minimization of in E 

For the numerical computations, we fixed D = 0.1, 7 = 1, oi = 02 = and 
/3i = /32 = 1 (Neumann boundary conditions). Besides, we assumed that 
Ui = 0.2, e = 0.3 and xq = 2/3. The integer n was set to 10 in the definition 
of E. 

Numerical computations were carried out for 100 functions ii}^m.E: 

n 

= J^/i' -i [{n - 2) (x - a)] ,k = i... 100, 

1=0 

whose components were randomly drawn from a uniform distribution in 
(-5,5). 

Minimizations of the functions G^^ were performed using MATLAB's® 
fminunc solver 0. This led to 100 functions fi^ in E, each one corresponding 
to a computed approximation for a minimizer of the function G^^.. In our 
numerical tests, we obtained values of G^^,{fil) in (5 • 10"'', 10~^), with an 
average of 5 • 10~^ and a standard deviation of 2 • 10~^. 

The values |l/^fc — M^||L2([o,i])/||Mfc||L2([o,i])) for A; = 1 . . . 100, are comprised 
between 5 • 10~^ and 0.16, with an average value of 0.04 and a standard 
deviation of 0.03. 

Fig. [11(a) depicts an example of function fi in E, together with a function 
fi* which was obtained by minimizing G^. 

4.3 Test of another criterion 

In this section, we illustrate that measurement ()2.9p alone cannot be used 
for reconstructing fi. Replacing G^ by: 

Hf,{fl) = \\u{-,Xo) - u{-,Xo)\\L\0,e), 

and setting H^{il) := ||n''(-,xo) — ^''^(•5 2;o)||L2(o,e)) we performed the same 
analysis as above, with the same samples /ifc € -E and the same parameters. 

The corresponding values of i/^(/x*) are comparable to those obtained 
in Section 14.21 Namely, these values are included in (2 • 10~*^, 10"^), with 

^ MATLAB's® fminunc medium-scale optimization algorithm uses a Quasi-Newton 
method with a mixed quadratic and cubic line search procedure. Our stopping criterion 
was based on the maximum number of evaluations of the function G^, which was set at 
2 • 10^ 
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0.2 0.4 0.6 0.8 1 

(a) /i (plain line) and ^* (dotted line) 
4 
2 


-2 
-4 
-6 

0.2 0.4 0.6 0.8 1 

(b) jj. (plain line) and /i* (dotted line) 

Figure 1: (a) An example of function /i in E, together with a func- 
tion /i* which was obtained by minimizing G^. In this case — 
A**IIl2([o,i])/|1/^|1l2([o,i]) = 0-03 and G^{fi*) = 2 ■ 10~^. (b) The same func- 
tion together with the function fi* obtained by minimizing iJ^. Here, 
11/^ - ^1Il2([o,i])/IImIIl2([o,i]) = 0.47 and Hf^i^j*) = 2 • IQ-^. 
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average 2-10~ , and standard deviation 3-10~ . However, the corresponding 
values of the distance \\nk — /^^||L2([o,i])/ll/^fc||L2{[o,i]) larger than those 

obtained in Section [4.21 these values are comprised between 0.08 and 1.64, 
with an average of 0.56 and a standard deviation of 0.34. 

Using the same sample /i G £" as in Fig. [T](a), we present in Fig. [I](b) the 
approximation fi* obtained by minimizing if^j. In this case, the distance 
Hl*\\l^{[o,i])/\\I-^\\l^{[o,i]) is 18 times larger than - I-i*\\l^([o,i])/\\i-i\\l^([o,i])- 

5 Discussion 

Studying the reaction-diffusion problem [V^^^y) with a nonlinear term of the 
type u{n{x) — 7tt), we have proved in Section [2] that knowing u and its 
first spatial derivative at a single point xq and for small times t G (0,e) is 
sufficient to completely determine the couple {u{t,x), n{x)). Additionally, if 
the second spatial derivative is also measured at xq for t S (0,e), the triplet 
x), //(x), 7) is also completely determined. 

These uniqueness results are mainly the consequences of Hopf 's Lemma 
and of an hypothesis on the set M of coefficients which fi{x) belongs to. 
This hypothesis implies that two coefficients in M can be equal only over a 
set having a finite number of connected components. 

The theoretical results of Section [2] suggest that the coefficients fi{x) and 
7 can be numerically determined using only measurements of the solution 
u of (V^^-y) and of its spatial derivatives at one point xq, and for t € (0,e). 
Indeed, the numerical computations of Section[4]show that, when 7 is known, 
the coefficient fi{x) can be estimated by minimizing a function G^. The 
function u being the solution of iVfi^^)-, we defined G^(/u) as the distance 
between {u, du/dx){-, xq) and {u,du/dx){-,XQ), in the L^(0,e) sense. 

The numerical computations presented in Section 14.21 were carried out 
on 100 samples of functions fik chosen in a finite-dimensional subspace of 
M. In each case, a good approximation of fi^ was obtained. The average 
relative L^-error between 1^1^ and /.t^ is 30 times smaller than the average 
relative L^-error between fi^ and the constant function /Ufc(xo). Thus, a 
measurement of u and of its first spatial derivative at a point xq (and for 
t E (0,e)) indirectly gives more information on the global shape of fj, than a 
direct measure of fi at xq. These good results, in spite of the computational 
error, indicate L^-stability of the coefficient n with respect to single-point 
measurements of the solution u of ("^^,7) and of its spatial derivative. 

Proposition 12.31 shows that the uniqueness result of Theorem 12.11 is not 
true without the assumption ()2.10p on the spatial derivatives. This suggests 
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that measurement (j2.9p alone cannot be used for reconstructing fi. In Sec- 
tion 14.31 working with the same samples /i^ as those discussed above, we 
obtained approximations /x* of /i^ by minimizing a new function H^, which 
measures the distance between u{-,xo) and u{-,xq). The average relative 
L^-error between Hk and fik* was 14 times larger than the average relative 
L^-error separating /i^ and /i^. This confirms the usefulness of the spatial 
derivative measurements for the reconstruction of /U. 
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Appendix A: maximum principle 

The following version of the parabolic maximum principle can be found in 
[271 Ch. 2] and [29l Ch. 3]. 

Theorem 5.1. Let u e Cf{{0,T] x (xi, X2)) n C([0, T] x [xi,X2]), for some 
T > and xi,X2 e M. Let c{t,x) < G Cj'^/2([0'^] ^ [xi,X2]), for some 
r?G(0,l]. 

Suppose that — D—^ — c(x)u > for t G (0,T] and x G (xi,X2). 

at ox'' 

(a) If u attains a minimum m* < at a point {t*,x*) G (0, T] x (xi,X2), 
then u{t,x) = m* on [0, t*] x [xi,X2]. 

(b) (Hopf's Lemma) If u attains a minimum m* < at a point {t*,xi) 
(resp. {t*,X2)), witht* > 0, then either ^{t* ,xi) >0 (resp. ^{t*,X2) <0) 
or u{t, x) = m* on [0, t*] x [xi , X2] ■ 

(c) If u > 0, the results (a) and (h) remain true without the assumption 
c{t,x) < 0. 

An immediate corollary of this theorem is: 

Corollary 5.2. The solution u{t, x) of (V^^-y) is strictly positive in (0, +00) x 
{a,b). 

Proof of Corollary [5721 Assume that it exists {t*,x*) G (0, +00) x (a, b) 
such that u{t*,x*) < 0. 
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Set w{t,x) = ne~'^*, for A > large enough such that 

c{t, x) := fi{x) — 7n — A < in [0, t*] x [a, b]. 

The function w satisfies: 

dw d'^w 

Since «;(0, x) = Ui[x) > in (a, 6) and w{t*,x*) < 0, the function w admits 
a minimum m* < in (0, t*] x [a, b]. From Theorem l5.1l (a). and since Ui ^ 0, 
this minimum is attained at a boundary point: it exits t' € (0,t*] such that 
w{t',a) = m* < or w{t',b) = m* < 0. Without loss of generality, we can 
assume in the sequel that w{t',a) = m* < 0. Prom Theorem 15.11 (b), we 
obtain ^{t',a) > 0. Using the boundary conditions in problem (V^^-y), we 
finally get: 

^ ^ dw , , , 
aim* = Pi —it',a) > 0. 
ox 

Using assumption (j2.5p . we get a contradiction. Thus u{t, x) > in (0, +oo) x 
{a,b). The conclusion then follows from Theorem 15.11 (c). □ 



Appendix B: numerical solutions of (V^^^) and (P/i ^) 

The equations {V^^^) and (Vfi^y) were solved using Comsol Multiphysics® 
time-dependent solver, using second order finite element method (FEM) 
with 960 elements. This solver uses a method of lines approach incorporat- 
ing variable order variable stepsize backward differentiation formulas. Non- 
linearities are treated using a Newton's method. The interested reader can 
get more information in Comsol Multiphysics® user's guide. 
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